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Abstract 

We show that the probability distributions P n {q,p\ y) := |(n|p,g; y)| 2 , which are obtained 
from squeezed states, obey an interesting partial differential equation, to which we give two 
intuitive interpretations: as a wave equation in one space dimension and also as a pseudo- 
diffusion equation. We also study the corresponding Wehrl entropies S n (y), and show that 
they have minima at zero squeezing, y = 0. 


1 Introduction 

This talk is based mainly on a work which was done in collaboration with Salomon Mizrahi 
from Brazil. 

Squeezed oscillator states are defined in terms of the bosonic creation and annihilation op- 
erators, a* := ^-(x - £), and a := ^-(x + 5;), as follows: 

l*;0 * lPi«;0 where z : = (« + *p)/^ > (*) 

and |0) is the ground state of the harmonic oscillator. Both 'D and S are unitary operators, T) 
creates the coherent state, and is defined by 

Q 

T>(q,p) := exp[za t — z m a] = exp[ipx - q—] , (2) 

and 5(f) is the squeezing operator: 

5(0 := exp[^(£a t2 - fa 2 )] , ( 3 ) 

where £ is a complex variable. For £ = 0, we recover the ordinary (unsqueezed ) coherent states. 

The squeezed states satisfy the completness relation, / |p> 9; £)(p, £| = 1 ^ f°r every £. 

Therefore, 

J Pn(q,p\0 = 1 , where P„(fl,j>;0 : = l(p,9;fl n )| 2 > ( 4 ) 

where |n) is the number state. If we interpret the real parameters q and p as the position and 
momentum variables, then (4) allows us to interpret the non-negative functions P n as proba- 
bility distributions in the (q,p)-phase plane, for every n and (. 
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In this talk, I shall consider these P n for real values of the squeezing parameter £, which will 
be denoted by y. In particular, I shall show that the P n (q,P',y) satisfy the interesting partial 
differential equations (9) and (12) , to which two intuitive interpretations can be given. Finally, 
I shall show that the Wehrl entropy S n (y) (14) of the P n must have their minima at zero 
squeezing, y = 0. 


2 Explicit Form of the Distributions P n 


The distribution P n (q,P',0 '■= |(«b>9;0l 2 gives the probability of funding n bosons (photons) 
in the squeezed states |g, y; £). It is a physically important quantity, and it has been calculated 
by different methods. The dependence of P n (q,p; () on n was studied by Schleich and Wheeler 
[2]. For £ = y, the P n is given by the following complicated expression [1,3,7]: 


P n (q,p\y) ■= l(p,?;y|n)| 2 


2y/7 

2" n!( 7 + 1) 


\H n (2,r];w)\ 2 exp 


_il±2£! 

. 1+7 J ’ 


n > 0, 


(5) 


where 



r, 1 '■= 


1-7 
1 + 7’ 


and 


w := 


q + HP 
7 + 1 ’ 


( 6 ) 


and where H n (2,r];w) are the generalized Hermite polynomials ( QHV ), which are defined in 
terms of the raising operatores R(a, 0; x) = ax — 0-^ , as follows [1]: 


H n (a, 0\ z) = R n (a, /?; z) • 1 = 


[n/2] 

E 


#— 0 


n! 

(n — 2s)!s! 



(7) 


These polynomials are equal to the standard Hermite polynomials for a = 2 and 0 — 1. In the 
limit, 0 — ► 0, these H„(x) becomes simple powers of x: P n (a, 0,z) = a n x n . Therefore, in the 
limit of zero squeezing, 7 — » 1, we have 7 — ► 0, so that the above Q'HV ’s become simple powers 
of w. Thus, for y — ♦ 0, equation (5) gives the following well-known Poisson distribution of the 
unsqueezed coherent states: 


P 

P n (q,p;0) = — exp 



n > 0, 


where p 2 := q 2 -(- p 2 , 


( 8 ) 


When discussing probability distributions, it is useful to think of the regions that are surrounded 
by the equipotential curves, P n (g,p; y) = const. ; I shall call these regions potential regions. 
Thus, the potential regions of the above Poisson distribution P„(y,p;0) are concentric circles in 
the (q,p)-plane. But for y ^ 0, these regions will have approximately elliptical shapes, whose 
the major axes lie along the p-axis for y < 0 and along the q-axis for y > 0. These regions 
become more elongated in one direction and narrower in the other, as |y| increases. 


3 The Partial Differential Equation for the P n 

Since the integral (4) of the distributions P n (<hP',y) over the whole (q,p)-space remains con- 
stant under squeezing, it is useful to think of the change of P n (<l,P\y) as functions of y as a 
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redistribution of probability densities in phase space, which maintains the positivity condi- 
tion P„(g,p;y) > o for all y. This redistribution of the P„(g,p; y) is governed by the following 
interesting and amazingly simple partied differential equation: 

^Pn(?,p;y(7)) = \ ^( ? ,p;y(7)) , where 7 := e 2 *' . (9) 

This equation was originally obtained [1] by straightforward but lengthy differentiation of the 
expression (5) , and by using the following property of the QHV [lj: 

,,«.)■ ( 10 ) 

However, we can now derive it by two other more general methods [5], as reported in the 
summary section. 

4 Interpretation as Wave and Pseudo-Diffusion Equa- 
tions 

I shall now present two possible intuitive interpretations of the above differential equation: 


(I) D’Alembert or Wave Equation: The following is a new interpretation, which was 
not discussed in [1]: For a fixed squeezing parameter y, equation (9) looks like the wave 
equation for one space dimension q, if we think of the p variable in (9) as the time variable 

t: 


( d 2 1 d 2 \ 


1 d 

where p(q,t;y) = --^P n (q,t;y(y )) , 

(11) 

In this interpretation, the parameter 7 would then play the role of the speed of light c(n) in 
matter, which depends on the parameter y, similar to the dependence of c(n) on the index 
of refraction index n. If the P n are thought of as electromagnetic potentials $(q, t; y), then 
4^P n (tf,p;y( 7 )) will play the role of a time-dependent charge distributions —4np(q,t;y). 

(II) Pseudo-Diffusion Equation: By substituting ^ = 2e 2y ^ into (9) , we obtain a more 
symmetric differential equation for the P n : 




( 12 ) 


This equation is also new and permits a more pertinent intuitive understanding of the 
redistribution process of the P n , by comparing (12) with the diffusion equation in two 
dimensions [6]: 


jd l_ dl) 

W + Sp’J 


T(q,p;t ) , 


—T{q,p-t) = a 
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(13) 


where a is the diffusion coefficient. Equations (12) and (13) are similar, if we interpret 
the squeezing parameter y as the time variable. However, the two equations differ in two 
interesting aspects: 

(1) The sign in front of ^ in (12) is negative rather than positive. Such a “negative 
diffusion coefficient” leads to a infusion” rather than diffusion in the p-direction. Con- 
sequently, as y increases, the equi- probability curves, P n (<hP\y) = const move towards 
the origin along the p-axis, but away from the origin along the q-axis. Therefore, we ex- 
pect the probability regions to be concentric elongated “quasi ellipses” which axe extended 
along the p-axis for y — ► — oo. They become more and more circular as y approaches zero, 
and then stretch outwards dong the q-axis, as y -4 oo. For the above reasons, we shall 
call equations (9) and (12) “pseudo diffusion equation”. 

(2) The “diffusion coefficients” exp[2y]/2 and -exp[-2y]/2 and in front of and ^ in 

(12) depend on y. For y — ► -f oo, the term \c 2y ^sPn dominates the r.h.s. of (12) , whereas 
for y — » — oo, the second term dominates. This dependence on y can be given an interesting 
intuitive explanation: Let us consider the redistribution process when y is very large: In 
this case the probability densities P n (g,p;y) are extended in the q-direction and tightly 
squeezed or compressed in the p-axis, which makes it difficult to compress them further 
along the p-axis. For this reason the “infusion coefficient” becomes so small, namely 
ex exp[— 2yj. In contrast, the diffusion along the q-axis must become faster and faster, in 
order to diffuse all the incoming density flux from the other orthogonal p-direction, which 
is entering the cigar-shaped potential regions through their lengthy boundaries. 


5 The Wehrl Entropy for the P n 

A useful measure for the information content of the probability distributions P n (q,p\y ) is the 
Gibbs or Wehrl entropy [7], which is defined by 

S„(y) := -/ />„(?, P;v)lnP„(«,P;!/) ^ ■ (14) 

Because of the symmetry P n (q,p\-y ) = P n (p,q;y), the entropy (14) is even in y: S n (-y) = 
S n (y). Therefore, at y = 0 each 5 n (y) must have either a maximum or a minimum. We shall now 
argue that 5„(0) should correspond to a minumum: We assume that S„(y) does not oscillate 
as a function of y. Therefore, it is enough to argue that S„(y) grows with |y| for large values 
of |y| . For large positive y, equation (12) behaves essentially like a one-dimensional diffusion 
equation in the q-variable. But it is well-known that the solutions of diffusion equations lead 
to entropies which increase with time [6]. Therefore, the S n (y) must increase as y — ► oo. But 
since the S n (y ) are even in y, they must also grow as y — ► — oo. Hence, the S„(0) must lie at 
the bottom of the curves S(y) vs. y. 

Finally, we note that the von Neumann entropy S v jv(p) := — Tr(plnp) for the pure states 
p := |n)(n| must vanisch. In contrast, explicit calculations of the Wehrl entropies of the Poisson 
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distributions (8) shows that S n (0) > 1 for all n, in accordance with a conjecture by Wehrl [7], 
which was proved by Lieb [8]. 

To summarize this section: in contrast to diffusion equations, where the entropies of their 
solutions always increase with time, the entropies S n (y) for the solutions of the above pseudo- 
diffusion equation first decrease monotonically as y grows from — oc to zero, but then increase 
monotonically as y grows from zero to +oo. 


6 Summary and Outlook 

Two equivalent partial differential equations (9) and (12) were presented and then interpreted, 
as wave and as pseudo- diffusion equations. The probability densities P n (q,p\y ) (5) provide 
infinite number of their solutions. 

By the time of writing the present lecture notes, we succeeded in proving, by two general 
methods, that the expectation values (g,p; i\0\q,p\ £) of an arbitrary operator 0, satisfy a 
generalized version of the above partial differential equations, which also include rotations, i.e 
for the general squeezing f = re‘*. Interesting examples of 0 are the number operators N and 
N 2 ; their expectation values provide the simplest solutions of (9) and (12) . Also the projection 
operator |g, p; £)(g,p; £|, and consequently its Wigner function, satisfy these equations. 
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